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RESUMEN 

En este articulo se describe un modelo 2D sintetico desarrollado en OpenSees, de 
un muelle apoyado sobre pilotes, capaz de representar el comportamiento no 
lineal de la estructura ante cargas ciclicas. Se explica los nudos y elementos del 
pilote, asi como la section transversal que es modelada con fibras. La interaction 
suelo-estructura es considerada a traves de resortes discretos no lineales 
horizontales y verticales. Se propone un procedimiento para obtener la curvatura a 
lo largo del pilote, a partir de las curvas esfuerzo-deformacion del hormigon en las 
fibras extremas opuestas. La rotation plastica se obtiene mediante integration 
numerica de una region plastica definida en cada instante de carga, cuando se 
alcanza la curvatura de fluencia. Finalmente se calcula la longitud plastica en el 
cuerpo y la cabeza del pilote, donde esta ultima se obtiene por la contribution 
independiente de la flexion y la adherencia. 

PALABRAS CLAVE: Dano, curvatura de fluencia, longitud plastica, resortes no 
lineales, metodo de fibras. 


ABSTRACT 

In this article a simplified 2D model developed in OpenSees of a wharf supported 
on piles is described. This model is capable to represent the nonlinear behavior of 
the structure under cyclic loads. Nodes and elements of the pile are explained, and 
the cross section is modeled following the fiber method. Soil-structure interaction is 
considered by nonlinear discrete horizontal and vertical springs. A method is 
proposed for the curvature along the pile from stress-strain curves of the concrete 
at the opposite end fibers. Plastic rotation is obtained by numerical integration of a 
plastic region defined at each instant of loading, when the yield curvature is 
reached. Finally the plastic length in both body and head of the pile is calculated. In 
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the head of the pile, the total plastic length is obtained as a product of the 
independent contribution of bending and adherence. 


KEYWORDS: Damage, yield curvature, plastic length, nonlinear springs, fiber 
method. 


1. INTRODUCCION 

El Ecuador cuenta con cuatro autoridades portuarias: Manta, Guayaquil, 
Puerto Bolivar y Esmeraldas. Asi como tres puertos petroleros que son: Balao, El 
Salitral y La Libertad; ademas de once puertos privados (Camara Maritima del 
Ecuador, 2015). Estas terminales portuarias permiten el comercio nacional e 
internacional del pais. Sin embargo sus actividades podrian ser interrumpidas por 
los danos ocasionados durante un terremoto, debido a que estas se localizan en 
zonas de alto peligro sismico segun (NEC, 2015). En el caso de Esmeraldas, la 
ocurrencia de un evento sismico provocaria graves perdidas economicas para el 
Ecuador por su afectacion directa a la refineria y al puerto de embarque de 
petroleo en Balao (Calderon, 2009). 

En la actualidad se estan construyendo veintisiete puertos y facilidades 
pesqueras artesanales en las provincias de El Oro, Santa Elena, Esmeraldas y 
Manabi. El Telegrafo (2014), MAGAP (2015); asi mismo se proyecta el puerto de 
aguas profundas de Manta (APM, 2015). Por lo tanto es necesario un analisis y 
diseno apropiado de los elementos estructurales que conforman el puerto, para 
asegurar un adecuado comportamiento ante un sismo. 

La estructura tipica de estos puertos es una losa maciza soportada por 
pilotes. Por experiencia en otras zonas de alto peligro sismico, tales como la costa 
oeste de los Estados Unidos, Japon, Nueva Zelanda, Italia (Borg, 2007), se 
conoce que los elementos mas susceptibles a dano, y que requieren mas 
atencion, son los pilotes. En ellos, el dano puede localizarse en la cabeza junto a 
la conexion pilote-losa, y/o en el cuerpo del pilote (Goel, 2010). 

Existe ya la acumulacion de resultados experimentales sobre el 
comportamiento de los pilotes y la conexion pilote-losa a partir del ano 1988 (Joen, 
Park y Priestley, 1988), y en trabajos posteriores: Jellin (2008), Silva, Seible y 
Priestley (1997), Sritharan y Priestley (1998) hasta la fecha (Foltz, 2011). Ademas 
se han realizado modelos relativamente sofisticados usando elementos finitos 
(Shafieezadeh, 2011), y otros modelos simplificados como el propuesto en (Goel, 
2010). Este ultimo tipo de modelo tiene como ventaja su relativa precision y 
sencillez, lo que permite su difusion en el dia a dia del analisis y diseno de estos 
puertos. Sin embargo, es muy sensible a la variation de la longitud de las rotulas 
plasticas, por lo que actualmente se estan buscando nuevas alternativas para la 
description del dano en estructuras portuarias. 

En este trabajo se describe un modelo 2D de estructura portuaria que usa 
el metodo de las fibras, y que considera la interaccion no lineal suelo-estructura, 
para representar el comportamiento real de la estructura. Con este modelo se 
busca describir correctamente el dano en los pilotes expresado por una longitud 
plastica, que se obtiene por integration numerica de una region de curvatura 
plastica, definida segun los limites de fluencia de los materiales para cada zona 
del pilote (cuerpo y cabeza). Ademas en la cabeza del pilote se calcula la longitud 
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plastica producto de la contribution independiente de la flexion y la adherencia, 
debido a que el estudio de Jellin (2008), muestra el aporte de la adherencia en la 
conexion pilote-losa. 

El modelo 2D se desarrolla en OpenSees, un programa estructural de 
codigo abierto, que en comparacion a otros paquetes comerciales tiene como 
ventajas principales su versatilidad y estrecho contacto con el mundo academico, 
lo que entre otras cosas permite su continua actualization con tecnicas 
avanzadas, modificables por los usuarios para adaptarlas a sus propias 
exigencias. El lenguaje de programacion de OpenSees se denomina Tel (Tool 
command language, que usualmente se combina con el grupo de herramientas 
GUI (graphical user interface) Tk, por lo que se denomina Tcl/Tk. 

2. DESCRIPCION 

El modelo 2D consiste en una subestructura formada por tres pilotes 
unidos a una losa de cubierta, sometido a una carga ciclica lateral y una carga 
axial constante en cada pilote. Este modelo se basa en estudios realizados por 
Shafieezadeh y DesRoches (2010), que hacen una simplification del modelo 3D 
del muelle tipo, constituido por una losa maciza soportada por pilotes, a un modelo 
2D. A su vez en Caiza (2013), el modelo 2D con todos los pilotes de un eje del 
muelle, se simplifica a un modelo 2D con los tres pilotes mas cercanos a la costa. 
Debido a que Kawamata (2009) y Yang et al. (2012), demuestran que en un 
muelle el dano mayor ocurre en los dos primeros pilotes. 

Caiza (2013) indica que los dos pilotes mas cercanos a la costa tienen una 
luz libre mas corta, mayor rigidez y absorben mas fuerza sfsmica, con lo cual son 
mas susceptibles a sufrir danos, por ello se los denomina pilotes “sfsmicos”. Caiza 
(2013), Kawamata (2009) y Yang et al. (2012). Los pilotes mas alejados de la 
costa, son mas flexibles y absorben menos carga sfsmica, por lo tanto sufren 
danos menores y se les denomina pilotes “no sfsmicos” (figura. 1). 


Altura sobre la 



Figura 1. Esquema de la estructura portuaria para el modelo 2D. 


358 


Abel Zambrano, Pablo Caiza 


En la figura 2 se presenta una vista en planta de una parte de la estructura 
portuaria, donde se aprecia la losa de cubierta apoyada sobre los pilotes. Se 
observa el eje 2 de la estructura, en el cual se indica el espesor y el ancho 
colaborante de losa para cada vano. En los primeros modelos, la losa se discretiza 
en elementos tipo viga con propiedades geometricas segun los diferentes anchos 
de losa entre los ejes A y B, y B y C (linea azul), pero en la modelacion final se 
considera un ancho de losa constante de 3m entre los ejes A y B, y de 6m entre 
los ejes B y C (linea verde), porque debido a la gran rigidez de los elementos de la 
losa, los resultados fueron similares. 



Linea del area cooperante 

Linea de simplificacion para el 
modelo de la losa 

Figura 2. Vista en planta de la losa apoyada sobre los pilotes. 


3. GEOMETRIA Y MALLA DE LOS PILOTES 

El modelo 2D de estructura portuaria desarrollado en OpenSees, consiste 
en un archivo principal llamado “NLcontrol.tcl” compuesto por seis rutinas 
adicionales que seran descritas mas adelante. En este numeral se explica los 
datos iniciales que requiere el archivo “NLcontrol.tcl” para generar la geometrfa y 
malla de los pilotes. 

En la figura 1 se aprecia que los pilotes poseen una altura o longitud libre 
sobre el nivel del suelo de 1 m para los pilotes 1 y 2, y de 5.5m para el pilote 3, a la 
cual se denomina Li. Ademas los pilotes poseen una longitud embebida bajo la 
superficie del suelo denominada L2. Esta longitud esta dividida en varias capas 
(figura 1) y los espesores se los identifica de la siguiente manera: L3 es el espesor 
del relleno de roca, L4 es el de arena marina, y Ls es el de sedimento arcilloso. 

El sistema de referenda x-z de la estructura esta situado en el pie del 
pilote 1 de acuerdo a la figura 1 , de tal manera que la coordenada en x del pilote 1 
esta en x=0, la del pilote 2 en x=4.0m y la del pilote 3 en x=10.8m. El diametro del 
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pilote es D=0.63m y el pilote se divide en 84 elementos de 0.25m de longitud. La 
numeration de los elementos empieza de abajo hacia arriba, es decir desde el pie 
del pilote hasta la cabeza. 


4. PROPIEDADES DEL SUELO 

Cada pilote posee una longitud embebida bajo la superficie del suelo 
L 2 = 20 m para los pilotes 1 y 2, y de L 2 = 15 . 5 m para el pilote 3. Esta longitud L 2 esta 
formada por tres capas de suelo mostradas en la figura 1 . Las propiedades de la 
capa de relleno de roca y de arena marina, corresponden a materiales reales del 
Puerto de Los Angeles (Caiza, 2013), y se muestran en la Tabla 1. 

Tabla 1. Propiedades de los materiales granulares. 


Material 

Peso especifico 
Y (kN/m 3 ) 

Angulo de friction interna 
<p (grados) 

Relleno de roca 

18.5 

37 

Arena marina 

10 

38 


Con los datos de los materiales (Tabla 1), se utiliza la rutina 
“get_pyParam.td” desarrollada por McGann (2010) para determinar la resistencia 
ultima del suelo p u it y el desplazamiento yso (cuando p es 50% de p u it), que son 
parametros que definen las denominadas curvas p-y. En API (1987), Reese y Van 
Impe (2001), se encuentra informacion detallada acerca de las curvas p-y, y de los 
parametros del suelo p u it, yso. La rutina “get_pyParam.td” requiere los siguientes 
datos: 


• puSwitch. Especifica la forma de calculo de p u it y yso. Puede ser segun 
las recomendaciones de API (1987) 0 de Brinch Hansen (1961). Se utilizo 
la primera. 

• kSwitch. Define la variation con la profundidad del modulo de reaction 
del suelo. La option 1 considera una variation lineal (API, 1987), y la 
option 2 considera una variation parabolica modificada de API (Boulanger 
et al., 2003). Se escogio la option 1. 

• gwtSwitch. Considera la presencia del nivel freatico para el calculo del 
modulo de reaction del suelo (option 2), y si no existe nivel freatico 
(option 1). Se utilizo la option 2. 

La rutina “get_pyParam.td” se usa para calcular los parametros p u it y yso 
en las dos primeras capas de suelo del modelo 2D, por tratarse de materiales 
granulares; mientras que para la tercera capa (sedimento de arcilla) se trabaja 
directamente con los siguientes datos: p u it = 2270kN/m 3 y yso = 0.01 m. 


4.1 Resortes para representar el suelo 

La interaction suelo-estructura es uno de los aspectos mas importantes 
del modelo, para lo cual se considera que los elementos del pilote, discretos 
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verticalmente, estan apoyados sobre una cimentacion no lineal tipo Winkler. Para 
representar el comportamiento no lineal, el suelo se modela con un resorte 
horizontal p-y y un resorte vertical f-za lo largo de la longitud embebida del pilote, 
excepto en el pie del pilote donde se utiliza un resorte q-z en reemplazo del 
resorte t-z. 

En la figura 3 se aprecia el modelamiento de resortes empleado en el 
modelo 2D de estructura portuaria, implementado por McGann (2010) para la 
libreria de ejemplos de OpenSees. McGann desarrollo tres rutinas para calcular 
los parametros del suelo y utilizar los elementos de OpenSees “PySimplel”, 
“TzSimplel” y “QzSimplel” (Boulanger, 2003) , con la finalidad de modelar el 
resorte p-y, t-zy Q-z, respectivamente. Un dato importante para el resorte Q-zes 
el modulo de code en el pie del pilote, que tiene un valor de G=150000kPa. 

Los parametros del suelo para los resortes t-z y q-z se calculan con las 
rutinas “get_tzParam.tcl” y “get_qzParam.td”, siguiendo un procedimiento similar al 
descrito para la rutina “get_pyParam.td”. 


Altura sobre el 
nivel del suelo 



Resorte p-y 


Nivel del suelo 


Longitud 
embebida 
del pilote 


Nudo del 
pilote '*5 


L2 


Elemento 
del pilote 

Iguales 
grados de 
libertad 



Resorte t-z 

Nudo de resorte fijo 
Nudo de resorte transitorio 

Resorte Q-z 


Figura 3. Resortes usados en el modelo para representar el suelo. Fuente: 

(McGann, 2010) 


4.2 Elementos de longitud cero para los resortes 

En la figura 3 se aprecia que en cada nudo del pilote existe un resorte 
horizontal y uno vertical, donde el resorte horizontal se ubica en el nudo de resorte 
transitorio y el resorte vertical en el nudo de resorte fijo. Tanto el resorte horizontal 
como el resorte vertical incorporan diferentes comportamientos en un mismo nudo, 
razon por la cual se utiliza el elemento denominado de “longitud cero” para 
representar los dos tipos de comportamiento (horizontal y vertical) de una zona del 
pilote, al mismo tiempo que es independiente de la longitud. Entonces el elemento 
de longitud cero esta definido por dos nudos que tienen la misma posicion pero 
con diferentes materiales, donde en este caso el material es un tipo de resorte. En 
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el pie del pilote el elemento de longitud cero creado esta formado por los resortes 
“PySimplel” y “QzSimplel”, y en el resto de la longitud embebida del pilote esta 
formado por los resortes “PySimplel” y “TzSimplel”. 

El elemento de longitud cero conecta los nudos transitorios con los nudos 
fijos, donde los nudos transitorios estan libres en desplazamiento horizontal y giro, 
mientras los nudos fijos estan completamente restringidos. Despues a los nudos 
transitorios se les da los mismos grados de libertad que a los nudos del pilote, es 
decir desplazamiento horizontal y vertical, y giro. Para los pilotes 1 y 2 se crean 80 
elementos de longitud de cero que representan el suelo con resortes, y para el 
pilote 3 solo 62 elementos, porque la longitud embebida del pilote L 2 es 20m para 
los primeros pilotes y 15.5m para el ultimo. 

5. NUDOS Y ELEMENTOS DEL PILOTE 

5.1 Nudos del pilote 

Los nudos del pilote tienen libertad en desplazamiento horizontal, desplazamiento 
vertical y giro, y son 85 nudos en cada pilote, ya que existen 84 elementos. En la 
figura 4 se observa la numeration de los nudos y elementos del pilote, los nudos 
de los resortes, asf como el elemento rfgido que conecta el pilote con la losa y los 
elementos discretos de la misma. 


Elemento discreto 
de la losa tipo viga 



O Nudo del pilote # Nudo de resorte rotacional en la conexion 

• Nudo de resorte transitorio © Nudo en la losa 

# Nudo de resorte fijo 

Nota: Esta numeracion corresponde al pilote 1 . Para los demas pilotes la 
numeracion continua de manera similar. 


Figura 4. Tipos de elementos y nudos que intervienen en el modelo de los pilotes. 
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5.2 Elementos del pilote 

Se utiliza el modelo de fibras desarrollado por Caiza (2013), que incluye el 
elemento tipo viga-columna no lineal basado en fuerzas de OpenSees 
(“nonlinearBeamColumn”) para modelar el pilote, porque se ajusta al 
comportamiento real observado en los resultados experimentales. Este elemento 
es capaz de representar el comportamiento no lineal a flexion y carga axial. Los 
analisis estadfsticos de resultados y errores de los modelos empleados se 
describen en Caiza, Andrawes y LaFave (2010) y Caiza, Andrawes y LaFave 
( 2012 ). 


La integracion a lo largo del elemento no lineal de OpenSees se basa en la 
cuadratura de Gauss-Lobatto (Neuenhofer y Filippou, 1997), con un punto de 
integracion en cada extremo del elemento. Puntos adicionales no son necesarios 
debido al caracter lineal de la grafica de momentos. 

5.3 Resorte rotacional 

En la conexion pilote-losa se crea un resorte rotacional con un elemento 
de longitud cero que conecta el nudo en la cabeza del pilote con un nudo en la 
misma ubicacion usando una seccion tipo fibra que simule el deslizamiento de las 
varillas de acero con respecto al hormigon. Este resorte rotacional permite 
considerar la contribution de adherencia en la conexion pilote-losa. La conexion 
corresponde al elemento 85 en cada pilote de acuerdo a la figura 4. 

La generation de los nudos y elementos del pilote, los elementos de 
longitud cero para los resortes y el resorte rotacional en la conexion para cada uno 
de los pilotes, se lo realiza con la rutina denomina “NLl.tcl”. 

5.4 Seccion tipo fibra para los elementos del pilote 

La seccion transversal del pilote se divide en varias fibras en el sentido 
radial y angular (figura 5), y en cada una se obtienen los esfuerzos y las 
deformaciones. Esta discretization permite ademas, modelar distintos materiales 
en una misma seccion como son: el hormigon no confinado de recubrimiento, el 
hormigon confinado encerrado por los estribos, el acero de preesfuerzo y las 
barras de acero de refuerzo. 



Hormigon no 
confinado 
Hormigon 
confinado 

• Barra de acero 

* Acero de preesfuerzo 



ro = radio de la seccion (hormigon no confinado) 
rc = radio del hormigon confinado 
rci = radio a los cables de acero de preesfuerzo 
rc2 = radio a las barras de acero de refuerzo 


Figura 5. Seccion tipo fibra del pilote 
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Cada modelo se crea utilizando la librerfa de materiales uniaxiales 
disponibles en OpenSees, los cuales proporcionan information detallada del 
comportamiento no lineal de los materiales en comparacion con las normas de 
construction (Caiza et al., 2012). 

Se emplean diferentes secciones de fibras a lo largo de la longitud del 
pilote, para tener un mejor modelamiento del comportamiento “real”. En el caso de 
la section para la conexion pilote-losa, interviene las barras de acero de refuerzo, 
mientras que para el suelo no actua y es necesario en cambio considerar el acero 
de preesfuerzo. 

Las secciones tipo fibra del pilote se crean con la rutina “NLmaterial.tcl”, 
donde se ingresa toda la information necesaria en cuanto a dimensiones, 
divisiones de la fibra, recubrimientos y propiedades de los materiales. 

5.5 Propiedades de los materiales de la section transversal 

El pilote tiene un diametro de 63cm y es de hormigon preesforzado con las 
siguientes caracteristicas: 

• Resistencia a la compresion del hormigon no confinado f c =420kg/cm 2 . 

• Resistencia a la compresion del hormigon confinado f C c=620kg/cm 2 . 

• Resistencia a la fluencia de las barras de acero de refuerzo 

f y =5000kg/cm 2 . 

• Area de las barras de acero de refuerzo A s =8cm 2 , area total = 64cm 2 . 

• Resistencia a la fluencia del acero de preesfuerzo f P =1 7000kg/cm 2 . 

• Area de los cables de preesfuerzo A sp =1cm 2 , area total = 22cm 2 . 

Los valores de resistencia de los materiales corresponden a la resistencia 
esperada, denominada “real”. De esta manera para el hormigon la resistencia 
esperada es fce=1.3f c , y para el acero f ye =1.1fy, (Priestley, Calvi y Kowalsky, 
2007). 

5.6 Elementos rigidos entre la cabeza del pilote y la losa 

En la figura 4 se aprecia un elemento rfgido que conecta la cabeza del 
pilote con el eje de la losa, el cual se modela usando un elemento elastico tipo 
viga-columna de OpenSees. El comportamiento elastico asumido para este 
elemento, se debe a su gran rigidez. Esta es de un orden de magnitud mayor a los 
elementos estructurales adyacentes. La losa de cubierta a la cual se unen los 
pilotes se discretiza en 43 elementos de 0.25m y se modela usando igualmente el 
elemento elastico tipo viga-columna. Se asume un comportamiento elastico 
basado en estudios analiticos y experimentales previos (Shafieezadeh y 
DesRoches, 2010). 

La creation de los elementos rigidos que conectan la cabeza de cada 
pilote con la losa de cubierta y los elementos discretos de la losa, se realiza en el 
archivo principal “NLcontrol.td”. 




364 


Abel Zambrano, Pablo Caiza 


6. CARGA 

El archivo “NL2.txt” corresponde a la rutina utilizada por el programa 
principal “NLcontrol.td” para generar la carga axial y la carga ciclica que se aplica 
a los pilotes. Esta ultima consiste en pequenos pasos de desplazamientos de 
carga y descarga con incrementos de 0.001m en unos ciclos y con 0.0001m en 
otros ciclos, hasta completar desplazamientos picos de 0.013m (ciclo 1) y 0.195m 
(ciclo 8). Vease tabla 2. 

Tabla 2. Desplazamiento pico de cada ciclo de carga 


Ciclo 

Desplazamiento pico (m) 

1 

0.013 

2 

0.026 

3 

0.039 

4 

0.052 

5 

0.065 

6 

0.098 

7 

0.130 

8 

0.195 


Para la resolution del sistema de ecuaciones se usa el metodo iterativo de 
Newton clasico (OpenSees, 2013). Sin embargo, si este no converge, se prueba 
en primer lugar Newton con tangente inicial (Newton with initial tangent) y, 
finalmente, Newton con busqueda lineal (Newton with line search). 

7. SECCION CRITICA EN EL PILOTE 

En los pilotes denominados “sismicos” se podrian formar rotulas plasticas 
en la cabeza y en el cuerpo del pilote (Goel, 2010), por ello se determina la 
seccion critica del pilote, para verificar la posible formation de la rotula plastica en 
esta zona. La seccion critica se localiza a una profundidad de -2.5m (elemento 
70), de acuerdo al diagrama de momentos en el pilote para cada ciclo de carga 
(figura 6). 


Diagrama de momentos en el pilote en cada ciclo 



Figura 6. Diagrama de momento maximo a lo largo del pilote en cada ciclo 
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8. CALCULO DE LA LONGITUD DE LA ROTULA PLASTICA 

En este numeral se explica el calculo de la longitud plastica del pilote 
basado en la determination de una zona del pilote donde la curvatura sea mayor a 
la curvatura de fluencia. La curvatura se la obtiene mediante las curvas esfuerzo- 
deformacion en la cabeza y en el cuerpo del pilote. 

8.1 Curvatura 

La curvatura es la relation entre la deformation unitaria del material en la 
fibra extrema y su distancia al eje neutro. En la figura 7 se muestra la distribution 
de deformaciones en la fibra positiva y negativa de la section del pilote. 


£c 



a) Section del pilote b) Deformaciones unitarias c) Deformaciones unitarias 

de la fibra positiva de la fibra negativa 

£c = deformation unitaria del hormigon en la fibra positiva 
£s = deformation unitaria del acero en la fibra positiva 
scn = deformation unitaria del hormigon en la fibra negativa 
8sn = deformation unitaria del acero en la fibra negativa 

Figura 7. Distribution de deformaciones en la fibra positiva y negativa de 
la section del pilote 

A partir de la figura 7b se obtiene la ecuacion de la curvatura: 


La ecuacion (1) requiere conocer la profundidad del eje neutro c, asi como 
la deformation del hormigon e c y del acero e s . Sin embargo en este trabajo la 
curvatura se obtiene unicamente a traves de las deformaciones del hormigon en la 
fibra positiva y negativa, ya que el OpenSees proporciona information de las 
curvas esfuerzo-deformacion en cualquier sitio. 


Donde: e c es la deformation unitaria del hormigon en la fibra positiva; e C N 
es la deformation unitaria del hormigon en la fibra negativa; D es el diametro del 
pilote. 
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Usando la ecuacion (2), se determina la curvatura a lo largo del pilote, en 
la cual se analiza dos secciones en particular: la seccion en la cabeza del pilote y 
la seccion critica en el cuerpo del pilote. 


8.2 Fluencia en la cabeza del pilote 

Se determina el instante en el que se produce la fluencia de la seccion en 
la cabeza del pilote, que puede darse cuando la deformacion unitaria del hormigon 
alcanza una deformacion de 0.002 o cuando el esfuerzo de las barras de acero 
alcanza el esfuerzo de fluencia f y , lo que ocurra primero. Una vez determinado ese 
instante, se obtiene la curvatura de fluencia en la cabeza del pilote para la fibra 
positiva (fry , cp y negativa <p y ,cN. En los terminos <j) y ,cp y <t> y ,cN, la letra “C” se utiliza 
para indicar que la fluencia es en la cabeza del pilote; la denominacion “P” o “N" 
indica el sentido de la fibra, positiva o negativa. 

8.3 Fluencia en el cuerpo del pilote 

Se determina el instante en el que se produce la fluencia de la seccion 
critica en el cuerpo del pilote (elemento 70), que puede darse cuando la 
deformacion unitaria del hormigon alcanza una deformacion de 0.002 o cuando el 
esfuerzo del acero de preesfuerzo alcanza el esfuerzo de fluencia f P , lo que ocurra 
primero. Una vez determinado ese instante, se calcula la curvatura de fluencia en 
el cuerpo del pilote para la fibra positiva <p y ,BP y negativa <p y ,BN. La denominacion “B” 
viene del termino ingles “body” que significa cuerpo, porque la letra “C” ya fue 
utilizada para la cabeza del pilote; la denominacion “P” o “N” indica el sentido de la 
fibra, positiva o negativa. 


8.4 Longitud plastica en el cuerpo del pilote 

La longitud plastica L p se calcula con la ecuacion (3): 



La ecuacion (3), tiene implicita la idealization de la curvatura propuesta 
por Park y Paulay (1983), quienes consideran que la curvatura plastica <p P es 
constante en una zona plastica de longitud L p . Ademas en la ecuacion (3), 0 P es la 
curvatura maxima de la region plastica y 0 p es la rotacion plastica obtenida por 
integracion numerica de la curvatura. 

Para cada instante donde la curvatura en la seccion critica del cuerpo del 
pilote (elemento 70) sea mayor a la curvatura de fluencia, se debe conocer 
cuantos elementos hacia abajo y hacia arriba del elemento 70 superan la 
curvatura de fluencia. Es decir se delimita una region de la curvatura del pilote que 
supere la curvatura de fluencia, luego con esa region plastica a traves de un 
proceso de integracion numerica se determina el area bajo la curva que viene a 
ser la rotacion plastica. 

En este artfculo se utiliza el metodo del trapecio para el calculo de la 
rotacion plastica y se plantea la siguiente ecuacion: 
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6 p - Area - + 4>p(i+i) + 0p(i+ 2 ) -I ^ 0p(n- 1 ) + * dE ^ 


4>p(i ) ~ 0(t) _ <Ay 
i = 1,2, ...,n 

Donde: <p y es la curvatura de fluencia en el cuerpo del pilote; <p(i) es la 
curvatura total en el cuerpo del pilote que alcanza la curvatura de fluencia; <p P (i) es 
la curvatura plastica en el cuerpo del pilote; dE es el diferencial de integration, que 
viene a ser el tamano de los elementos del pilote 

En la figura 8 se muestra un esquema del calculo de la rotation plastica a 
partir de la curvatura, en la que intervienen las variables descritas en la ecuacion 
(4). Se aprecia un area sombreada desde la section /' hasta la section n que 
corresponde a la rotula plastica. Ademas se observa que la longitud plastica del 
pilote L p , es la altura del rectangulo equivalente que tiene la misma area que la 
distribution real de la curvatura (Park y Paulay, 1983). 



Figura 8. Esquema del calculo de la rotation plastica en el cuerpo del 

pilote 

En la figura 9 se resume en un cuadro de barras la longitud plastica en el 
cuerpo del pilote para cada ciclo, donde se aprecia que la longitud plastica maxima 
es de 0.87m y se presenta en el ciclo 6 para un desplazamiento de 0.098m (Ver 
tabla 2). La primera rotula plastica se desarrolla en el ciclo 4 para un 
desplazamiento de 0.052m. 
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1 2 3 4 5 6 7 8 

Ciclos 


Figura 9. Longitud plastica en el cuerpo del pilote para cada ciclo 


8.5 Longitud plastica por flexion en la cabeza del pilote 

La longitud plastica por flexion en la cabeza del pilote se determina usando 
la ecuacion (3), que requiere como datos: la rotacion plastica 9 P y la curvatura 
plastica cp v . Utilizando la ecuacion (4), se obtiene la rotacion plastica, mediante 
integration numerica de la region de curvatura plastica. Esta region plastica es 
definida para cada instante en el cual, la curvatura (p (i) en la cabeza del pilote es 
mayor a la curvatura de fluencia 4> y (figura 10). Por lo tanto esta region empieza en 
la cabeza del pilote (seccion 84) hasta una seccion debajo de la cabeza (seccion 
84-n) que cumpla la condition de que $ (i ) > $ y . Notese que en la ecuacion (4), la 
curvatura plastica es: 4> p ^ = $ (i) - 4> y > 0. El valor de 4> p requerido en la ecuacion 
(3) es la curvatura plastica maxima, que siempre se localiza en la cabeza del pilote 
(seccion 84). 


Los resultados de la longitud plastica por flexion en la cabeza del pilote se 
muestran en la figura 1 1 . 



Figura 10. Esquema del calculo de la rotacion plastica en la cabeza del 

pilote 
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1 2 3 4 5 6 7 8 

Ciclos 


Figura 11 . Longitud plastica en la cabeza del pilote para cada ciclo 


8.6 Longitud plastica por adherencia en la cabeza del pilote 

La longitud plastica por adherencia se determina usando la ecuacion (3), al 
igual que para la flexion. Pero en este caso la rotacion plastica G p , se obtiene 
directamente de los archivos de salida de OpenSees en el nudo 247 mostrado en 
la figura 12. 


Losa 

Resorte rotacional 
0.25m 




Nudo 247 

Nudo 246 
Pilote 


Figura 12 . Nudos en la cabeza del pilote 1 

La rotacion por adherencia B a en la cabeza del pilote viene dada por la 
rotacion del resorte rotacional d 2 47 , es decir 0 a = d 2 47 . Obtenida la rotacion e a , se 
define la rotacion de fluencia en la cabeza del pilote 9 y de acuerdo al numeral 8.2. 
La longitud plastica por adherencia se calcula usando la ecuacion (5). 


En la figura 13 se muestran los resultados obtenidos de la longitud 
plastica por flexion y por adherencia en la cabeza del pilote, donde se aprecia que 
la contribution por adherencia es mayor a la flexion. 
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Figura 13. Longitud plastica por flexion y por adherencia en la cabeza del pilote 


8.7 Longitud plastica propuesta por Priestley et al. 1996 

La ecuacion (6), propuesta por Priestley et al. (1996), permite determinar 
la longitud de la rotula plastica en la cabeza del pilote. 

L p = 0.08 L t + 0.022 f y d bl (6) 

Donde: L t es la distancia desde la rotula plastica al punto de inflexion; f y 
es el esfuerzo de fluencia de las barras de acero en MPa; d m es el diametro de las 
barras de acero (refuerzo longitudinal) en metros. 

El primer termino de la ecuacion (6), representa la propagation de la 
plasticidad debido a la longitud del miembro (flexion), y el segundo termino 
corresponde a la longitud adicional asociada a la penetration de deformation en la 
losa (adherencia). 

La longitud L t es aproximadamente igual a 1 ,32m. Esta longitud se obtiene 
localizando el punto de inflexion en el diagrama de momentos a lo largo del pilote 
(figura 6). 

El esfuerzo de fluencia del acero de refuerzo es f y = 5000 kg/cm 2 y el 
diametro del acero de refuerzo es de d bl = 0.032 m, de acuerdo a las propiedades 
descritas en el numeral 5.5. 

Reemplazando los datos en la ecuacion (6), se obtiene que el primer 
termino correspondiente a la flexion es igual a 0.11m, y el segundo termino 
correspondiente a la adherencia es igual a 0.34m. Estos valores se resumen en la 
Tabla 3 y se comparan con los resultados obtenidos en el modelo 2D. 
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Tabla 3. Longitud plastica en la cabeza del pilote 


Longitud plastica 

Ecuacion de Priestley (m) 

Modelo 2D (m) 

Por flexion 

0.11 

0.12 

Por adherencia 

0.35 

0.66 

Total 

0.46 

0.78 


9. CONCLUSIONES 

El modelo 2D de estructura portuaria descrito en este artfculo, es una 
alternativa novedosa para la modelacion de un muelle, porque trata de representar 
de una manera mas cercana a la realidad su comportamiento estructural 
considerando el efecto de la interaction suelo-estructura y las propiedades de los 
materiales de la section transversal a traves del metodo de fibras. Asf mismo este 
modelo reduce el tiempo de computation porque analiza unicamente los tres 
primeros pilotes cercanos a la costa, puesto que son los mas propensos a sufrir 
dano. 


Se confirma que en un pilote cercano a la costa, denominado en este 
trabajo como pilote “sismico”, el dano puede suceder en la cabeza y/o en el 
cuerpo del pilote. De hecho en el pilote 1 del modelo 2D, la fluencia en el cuerpo 
del pilote ocurre en el ciclo 4, mientras que la fluencia en la cabeza del pilote tuvo 
lugar en el ciclo 6. Por lo tanto la primera rotula plastica se forma en el cuerpo del 
pilote, alcanzando una longitud de 0.87m. 

En este trabajo, en la cabeza del pilote se considera no solo el dano 
causado por la flexion, sino por la adherencia, la cual es modelada con un resorte 
rotacional que representa el deslizamiento de las varillas con respecto al hormigon 
en la conexion pilote-losa. La longitud plastica en la cabeza del pilote se obtuvo 
mediante la contribution independiente de la flexion (0.12m) y la adherencia 
(0.66m), llegando a un valor total de 0.78m. 

La Tabla 3 muestra que la longitud plastica por adherencia es mayor que 
la de flexion, ya sea utilizando la ecuacion (6), (Priestley et al., 1996) o con el 
modelo 2D. La contribution de adherencia en la longitud plastica total es de 76% 
para la ecuacion (6), y de 85% para el modelo 2D. Esto indica que la contribution 
por flexion y por adherencia en ambos casos, guardan similitudes en proportion. 
Sin embargo los valores obtenidos con la ecuacion de Priestley et al. (1996) son 
menores que los del modelo 2D, principalmente la longitud plastica por 
adherencia. Esto demuestra que los metodos simplificados pueden subestimar la 
longitud de la rotula plastica obtenida con el modelo 2D. 
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